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We present arguments that the CKM CP-violation in the standard model may be 
sufficient for the generation of the baryon asymmetry, if the electroweak transi- 
tion in the early universe was of the cold, tachyonic, type after electroweak-scale 
inflation. A model implementing this is described which complies with the CMB 
constraints and which is falsifiable with the LHC. Numerical simulations of the 
tachyonic transition with an effective CP-bias indicate that the observed baryon 
asymmetry can be generated this way. 
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1. Introduction 

(D ! 

(— I , There are three families in the minimal Standard Model (SM), why? Such 

why-questions lie outside the scope of the model, but the thought 'to al- 
low for CP-violation and baryogenesis' comes up easily. What is currently 
known experimentally about CP violation^ can be attributed to the CKM 
' matrix being complex, which is only possible for three or more families. 

This provides a strong motivation for finding a scenario for baryogenesis 
that makes use of this fact. SM electroweak baryogenesis^, is believed by 
many to be impossible, because of (?) the weakness of CKM CP-violation, 
and (a) the smoothness of the finite-temperature electroweak transition'^. 

Of course, the SM assumes neutrinos to be massless, and a better frame- 
work is its extension that includes right-handed neutrino fields with renor- 
malizable Yukawa couplings and a Majorana mass matrix, which may be 
called the Extended Standard Model (ESM). Leptogenesis is an attractive 
possibility within this framework^. However, it makes use of physics at the 
scale of 10^° GeV and it seems worth putting effort into scenarios closer to 
what we know on the scale < 200 GeV. So we must face issues (i) and {ii). 

In the following a scenario^'^''''^'^ will be reviewed based on a cold 
('tachyonic') electroweak transition caused by the coupling of a gauge- 
singlet inflaton to the ESM. First, we shall argue that CP violation is much 
stronger at zero temperature than often stated. 
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2. Magnitude of CKM CP-violation 

The usual estimate for the magnitude of CP violation is^'^'^ 

^cp = J {ml-ml){ml-m'^){m^-ml){ml-ml){m'l-ml){ml-ml)/T^'^, 

(1) 

where m„, . . . ,mf, are the quark masses and^ 

J = \lm{VfgVh^V;,V:g)\ - (2.88 ± 0.33) x lO^^ (2) 

is the simplest rephasing-invariant combination of the CKM matrix V. For 
T = 100 GeV, of order of the finite- temperature electroweak phase transi- 
tion, (1) gives the discouraging value Scp ~ 10^^^. At zero temperature one 
would replace nij/T by the Yukawa coupling Xj = \/2mj/v, with v = 246 
GeV the vev of the Higgs field, giving 

5cP = J(A2-A2)(a2-A2)(a2-A2)(a2-A2)(a2-A2)(a2-A2) « (3) 

even smaller than (1). However, even with the Higgs field settled in its vev, 
the measured CP violating effects in accelerator experiments are at a much 
higher level than 10"^'^, e.g. the magnitude of the decay asymmetry e' in 
the K'^ — system is about^ 4 x 10^*^. This suggests that the above order 
of magnitude estimates of 5cp are misleading, at least at zero temperature. 

To make the discussion more concrete, consider the effective action ob- 
tained by integrating our the fcrmions, Tdi = — Tr(ln_D), with D the Dirac 
operator. In studies and discussions of electroweak baryogcnesis, CP viola- 
tion has been taken into account by approximating the CP asymmetry in 
Feff by a leading dimension-six term^^ 



S^^VMF-'-i^,.), (4) 



where -F^i/ is the SU(2) field strength tensor, F^^ its dual, and M is a mass 
depending on the scale of the problem. It could be the mass scale of an 
extension of the ESM. Within the ESM, at finite temperature, one could 
take M = T, but what to use when T = 0? The fermion masses are cx {(p), 
but (fi^ (fi / {(fi^ (fi) does not seem to make sense. The ordering of the terms in 
Feff according to increasing dimension is questionable in a zero-temperature 
transition in which (<p^(^) increases from zero to some finite value. 

Luckily, there exists a very detailed calculation of the zero-temperature 
effective action in a general chiral gauge model, based on a gauge-covariant 
derivative-expansion that is completely non-perturbative in the Higgs field, 
by L.L. Salcedo^2,i3_ have apphed his results to the ESM (without 
Majorana mass terms), and found that (4) is incorrect for this case: the 
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rephasing invariant J does not appear as a coefficient of tT{F^'^ F^^) times a 
function of the un-diffcrcntiatcd Higgs field The first CKM CP-violating 
contribution lies unfortunately still beyond the scope of Salcedo's results, 
but a typical term is expected to have the form, in unitary gauge — > W^, 

where a,b,c,d are SU{2) indices and n'^bcdi^^) ^ non-trivial function 
of the Yukawa couplings times the Higgs field, dj = Xj h. For e.g. 

a = b = c = d= l the above expression violates CP and is expected 
to contain the invariant J. The crucial point is now that n^^^^ is expected 
to be a homogeneous function of the dj. This is the case for the explic- 
itly calculated coefficient functions at fourth order of the gauge-covariant 
derivative-expansion (involving four Lorcntz indices), which arc homoge- 
neous of degree zero^^'^'^, which is why (4) cannot occur. For the higher 
order term we expect that, when we rescale the Yukawa couplings Xj — > sXj , 
the coefficient functions n{Xh) scale by some negative power of s (perhaps 
up to logarithms). This strongly suggests that we should not include the 
product of Yukawa couplings (3) in estimating (5cp, leaving only J, which 
is a factor of about 10^'' larger. For example, 

{XI - Xl){Xl - Af)(Af - Xl){Xl - Xl){Xl - XDiXl - A^) . , „ 
[Xl + Xl){Xl + A?)(A? + XDiXl + A2)(A2 + Xl){Xl + Xj) ' ' ' 

The argument applies only at zero temperature, since at finite T a new 
scale appears, e.g. the thermal QCD quark mass mth = gsT/V&. 

3. A falsifiable model of electroweak-scale inflation 

Here we take seriously a minimal phenomenological extension of the ESM 
with an additional gauge-singlet inflaton that couples only to the Higgs 
field^. Its inflaton- Higgs potential is constructed after 

V{a, 4>)=Vo- ^a^aP + -^a.a" - ^X^^a^cf^ + i/z^,/,^ + -^X^{4,^f, (5) 

where a is the inflaton field and (f)^ = iLp^ Lp. During inflation, a slow-rolls 
away from the origin. After inflation has ended it accelerates, and the 
effective Higgs mass 

/^eff = /^l - Aa^O-^ (6) 
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changes sign from positive to negative, and the tachyonic electroweak tran- 
sition starts. During the transition the baryon asymmetry is to be gener- 
ated through the electroweak anomaly and the CKM CP-bias. To avoid 
sphaleron washout of the asymmetry the reheating temperature Trh has 
to be small enough. Taking Vf^^ = 100 GeV, the Hubble rate is tiny, 
H « lO-i-* GeV, and T^h = [30K)/(7r^5*)]^/'' - 51 GeV (5, = 86.25 is the 
effective number of SM degrees of freedom below the W mass), at which 
temperature the sphaleron rate is negligible. 

As a phenomenological model it has to comply with what is cur- 
rently known. Firstly, we have the CMB constraints^^: the large-scale 
normalization |^k|^ = (3.8 ± 0.5) x 10~^, and the scalar spectral-index 
ns-1 = -0.03±0.G3. The latter is given by - 1 ~ -{2/N){p-l)/ {p-2), 
with A'^ the number of e-folds before the end inflation, and > 2 for slow- 
roll to end naturally. For electroweak-scale inflation N c± 24, which is much 
lower than usually mentioned « 60, and the resulting ng — 1 < —0.083 is 
too low. This is fixed^ as mentioned below. 

Secondly, in the minimum of the potential we have (approximating the 
cosmological constant by zero) 

^{va,V4,)=0, ^{va,v<j>) =0, V{v„,v^)=0, = 246 GeV, 
aa acp 

t^l - K4>vl = -\ml, ml = 2X^vl, (7) 

where would be the SM Higgs mass, were it not that Xacp leads to a 
considerable inflaton-Higgs mixing, to be addressed shortly. For definiteness 
we continue with V^^'^ = 100 GeV, = 100 GeV, = 200 GeV, for 
which the Higgs self-coupling = 0.33, and we expect a particle mass 
above the current lower bound of 114 GeV^. 

Thirdly, the transition has to be reasonably rapid to allow for successful 
baryogenesis, which leads'''^ to the constraint Xa^, > 0.01. Essentially this 
constraint now fixes p = 5 (the choices for Vq, iJ.<p and are not critical 
in this respect). The basic reason is that for p < 5 the potential goes down 
rather slowly, giving a large cr-vev Va- and consequently a too small A^^ (cf. 
(7)), whereas for p > 5 the result is a ridiculously large Xa,p- So we end up 
with non-renormalizable couplings, and to keep the powers of a as low as 
possible we have chosen q = 6. 

Returning to the minimum of the potential, the inflaton-Higgs mixing 
follows from the mass matrix 



d^V/d{a, 4))v„,v^ ^ TOi = 385 GeV, wa = 125 GeV, ^ = 0.44 , 
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where ^ is the mixing angle. The model thus predicts (only) two scalar 
particles with couplings to the rest of the SM determined by the mixing. 

The CMB scalar spectral-index can be accommodated by adding terms 
— ^a2(T^ + a^a^ — ^a^cr'^ to the potential that were artificially set to 
zero in (5). For example, still with as = we found that we could fit ns 
easily, see figure 1. Furthermore, in a one-loop study we verified that the 




0.06 0.08 0.12 0.14 0.06 0.08 0.12 0.14 



Figure 1. Dimension-full (left) and dimension-less (right) parameters as a function of 
X = "2nipj^j^^j^/Vb, for the central WMAP value ng ~ 1 = —0.03. 

potential can be kept flat enough so as not to upset slow-roll inflation; the 
non-renormalizable couplings lead to the conclusion that the model breaks 
down at scales around a TeV or so. 

The flatness of the potential is finely tuned, but we stress that there 
is nothing wrong with this in a phenomenological model. It is a mini- 
mal extension of the SM that incorporates inflation as well as the tachy- 
onic electroweak transition desired for baryogenesis, with signatures due to 
inflaton-Higgs mixing that allow it to be falsified by the LHC. 

4. Tachyonic way-out-of-equilibrium transition 

After inflation and before the transition, the universe is cold and approx- 
imately in the ground state corresponding to the symmetric phase of the 
SM. Then /i^ff changes sign and the tachyonic transition develops. The 
initial part of such transitions with given 'speed' dfi^g / dt has been studied 
analytically as well as numcrically^^'^^'^®'^^'^°. We used the infinitely-fast 
quench, 

Mcff = +M^ i < 
= -p?, t > 0, 
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in our simulations of the S'[/'(2)-Higgs model described by 

as an explorative approximation for studying baryogcncsis^. In the gaussian 
approximation the Fourier modes of a real component (j)]^ of the Higgs field 
can be solved in the form ake""^** + f3ke'"^i'*^ with = y^— /x^ -|- fc2_ 
The modes with k < fi arc unstable and grow exponentially cx el'^fc After 
some time jit ^ 1, the equal-time correlators of (pu and TTk = dt4>k are 
generically well approximated by the classical distribution^^ 



exp 



2 „ , ynk + l/^ + fik nk + l/2-nk 



(9) 



where = , (wfe 0k ± TTk), and nfe + l/2 + nfe ^ 1 and nj; + 1/2 — 
rife <C 1. So the typical grow large, while ^ 0; the distribution gets 
squeezed. In the classical approximation**, realizations of the distribution (9) 
are used as initial conditions for time evolution using the classical equations 
of motion. For a weakly coupled system such as the SU(2)-Higgs sector of 
the SM we don't have to wait till the occupation numbers have grown 
large, as the quantum and classical evolutions are identical in the gaussian 
approximation. Our 'just a half initial conditions start directly at t = 0, 
for which the occupation numbers Uk = fik = 0. In temporal gauge, Aq = 0, 
the gauge field A is initialized to zero, and the initial SjA is found from 
Gauss' law^. 

In wo studied the particle-number distribution of the Higgs and W. 
For example. Fig. 2 shows how the W-particle rife in Coulomb gauge change 
in time. The rapid rise of the low momentum modes end at about truH = 
7 and after tmn = 30 the distribution does not change much anymore. 
These rife differ initially very much from the corresponding ones in the 
unitary gauge, but by time tniH = 100 they are practically the same. The 
particle numbers of the longitudinal polarization are initially similar to 
those of the Higgs but later appear to be slower in adjusting to the same 
distribution. The distributions at tniH = 100 can be fitted to Bose-Einstein 
distributions, separately in the more important low momentum region and 
for the higher momentum tails, giving effective temperatures and (at that 
time still c,onsidera})l{^) chemical potentials. The overall lesson learned from 



^Scc "^'^ for a comparison in the 0(4) model with the quantum case using ^-derivable 
approximation. 
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Figure 2. Coulomb-gauge W-particle numbers vs. momentum, time in units of frij^ . 

these studies is that the transition is strongly out of equihbrium, and over 
after a few tens of roj^ . The finer details of thermalization will be different 
when the other d.o.f. from the quarks, gluons, etc. have been taken into 
account. 



5. Baryogenesis 

To estimate the asymmetry generated by the transition under influence of 
CP violation we simulated^ the S'C/(2)-Higgs model with the effective CP 
violating term (4) added to the lagrangian (8). Although it was mentioned 
in section 2 that this term does not apply to the ESM, it is still interesting 
to see how large its parameter (5cp has to be in order to obtain the observed 
asymmetry. For the simulation we relabelled 

^ K, k = l&TT^m^ K. (10) 



167r2Af2 

The baryon asymmetry is given by the anomaly equation 

B{t) - 3Ncs{t) - 3^ f dt' I d^x (tr (^F^,{x,t')F^'' {x,t')')) 



167r2 Jo 

where we assumed B{0) — (after inflation, just before the transition) and 
initial Chern-Simon number A^cs(O) = (choice of gauge). Another useful 
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observable is the winding number of the Higgs field, 

= 2^ I d^XEMmtridkUU^dlUU^d^UU^), 

At relatively low energy the gauge and Higgs fields are close to being pure- 
gauge and then iVcs(^) ~ N^{ip). 

We simulated with parameters such that mu/rnw = 1 and v^, and 
various values of k defined in (10). The initial conditions were 'just a half 
(cf. Sect. 4) and 'thermal'. The latter are simply Higgs- field realizations 
drawn from a free BE ensemble at low temperature T = 0.1m//, thermal 
noise that we used to get an impression of the sensitivity to the initial con- 
ditions. Figure 3 shows an example of how at first {(f^ip) rises exponentially 




m„t 



Figure 3. Averages: Higgs field 2{(p^ (p) /v^ , Chern-Simons number (A^cs) and Higgs 
winding number (Aw), for mu = "Iid, k = 16, 'just a half initial conditions. 

and then settles into a damped oscillation near the vev. The Chern-Simons 
number has a small bump near tmn = 7, which we can understand ana- 
lytically, after which it appears to resonate with {(f^f) and settle near the 
winding number. The latter is at first erratic but then appears to stabilize 
earlier than iVcs- 

Fig. 4 shows results for the Chern-Simons number-density at time 
tniH = 100. There is no evidence for a sensitive dependence on the initial 
ensemble, and apparently the dependence on the magnitude of CP viola- 
tion becomes non-linear for fc > 5 — 10. The dependence on the mass ratio 
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mH/rnw appears limited. This is in contrast to what we found^"'^ in the 
1+lD abelian-Higgs toy model, for which there is a sensitive dependence 
on m^f/mvy; it also showed nice linear behavior in the 1+lD analog k for 
small enough values. The dashed lines are linear fits through the origin, 




Figure 4. Chern-Simons number density vs. k. 



ignoring the data at fc = and 16. The straight full lines go through the 
two data points with lowest k. Including the data at fc = is meaningful 
because these were generated with the same sequence of pseudo-random 
numbers as for k = i {nin = and k — 1.5 {mn — \/2mw)- (Note 

that the k — data is shifted somewhat away from zero for mn = mw)- 
Using the dashed-hne fit for the case uih = V^mw — 114 GeV gives a 
baryon asymmetry 

— = (0.46 ±0.21) X IQ-'^nmL. 

To fit the observed UB/n-y ~ 6.5 x 10^^", requires k to be 

_ 1.4 X 10"'^ _ 2.2 X IQ-^ 
'^^ ^ ITeV^ 

The 1 TeV scale is reasonable if we interpret Cqp as being due to physics 
beyond the ESM. On the other hand, if we boldly assume that the results 
make sense for the ESM by taking the scale M = mw, then the required 
6cp turns out to be of the order of J in (2) : 

Sep ~ 0.7 X 10"^ 
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We consider this as very encouraging for further developing the scenario of 
cold electroweak baryogenesis. 
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